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Abstract 

We show that the variation of the ground state entanglement in hnear, higher spatial derivatives 
field theories at zero-temperature have signatures of phase transition. Around the critical point, 
when the dispersion relation changes from linear to non-linear, there is a fundamental change in 
the reduced density matrix leading to a change in the scaling of entanglement entropy. We suggest 
possible explanations involving both kinematical and dynamical effects. We discuss the implication 
QJ of our work for 2-D condensed matter systems, black-hole entropy and models of quantum gravity. 
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I. INTRODUCTION 



A composite quantum system can be prepared in the so-called Entangled (inseparable) 
states. These states precisely measure the non-local quantum correlations between the 
subsystems which can not be explained by classical physics. Entanglement is at the heart 
of some of the deepest puzzles of quantum mechanics and its most promising applications 
[1, 2]. 

In recent years, entanglement entropy is found to be playing crucial roles in understanding 
quantum behavior of macroscopic and microscopic systems. For the macroscopic systems 
like black-holes, entanglement entropy refers to the measure of the information loss (for an 
outside observer) due to the spatial separation between the degrees of freedom inside and 
outside the horizon. Interestingly, the entanglement entropy turns out to be proportional to 
the area of the horizon [2-5] and raises the possibility to interpret the Bekenstein-Hawking 
entropy as the entanglement entropy. 

At the microscopic level, it has been observed that the entanglement plays an important 
role in quantum phase transitions [6-8]. Typically at the critical point, where the transition 
takes place, long-range correlations in the ground state develop. Like the free energy for the 
ordinary phase transitions, entanglement is an useful quantifying measure to understand the 
long-range correlations in quantum phase transitions [9-13] . The simplest versions of phase- 
transitions assume the (scalar) order parameter to have non-linear self-interactions F{(f)) — 

a20^ -I- 040^ H [14]. One way to generalize the above model is to include inhomogeneities 

in the order-parameter of the form [15]: 

oo oo 

Fid^) = J2 «2n<^'" +Y.bm {V^4>f . (1) 

n=l m=l 

We show that a fundamental change in the ground state of a quantum system arises due to 
higher spatial derivative terms (especially, 63 7^ 0) even if the self-interactions are set to zero 
(a2n = 0). In particular, we show that for a real scalar field (0) whose Hamiltonian is [16-19]: 



H = ^ I (fx 



7r2(x) + f]6^(V^ 



m=l 



(2) 



the entanglement entropy, computed for the discretized Hamiltonian, shows a sudden dis- 
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continuity due to odd coupling coefficients {bmY- We explicitly show this up to m = 3. In 
the above Hamiltonian, we let bm — ±k~^^'^~^^ where k is the spatial frequency at which 
the dispersion relation changes from linear to non-linear. In the next section, we briefly 
discuss the field theory model, we consider, and the computational method to determine the 
entanglement entropy. 
II. THE FIELD THEORY MODEL 

Our model is the Hamiltonian (2) which corresponds to a real scalar field propagating in 
(3 + l)-dimensional flat space-time with higher spatial derivative terms at zero temperature. 
To extract the physics, we will consider the following dispersion relation: 

,;^=e+'k'+l^k\ (3) 

keeping the terms up to m = 3 in Hamiltonian (2), where e, r are (dimensionless) constants. 

It is harder to obtain an analytical expression for entanglement for free fields [3, 4]. In 
this work, we use the following semi-analytic approach: 

(a) Discretize the Hamiltonian (2) along the radial direction of a spherical lattice with lattice 
spacing a, such that r ^ Vi — ia; r^+i — ri — a: 

1 ^ 

^ = 5Z = ^ H ^^J-^l + '^3 Kij ■ (4) 

3 i,3 

Note that (/7j's, tt^'s and K^j^s are dimensionless [4] and the interactions are contained in 
the off-diagonal elements of the matrix K^j (see Appendix A) that are functions of P, £, N, e 
and T^. Here P = and P'^-^ is the (dimensionless) effective couphng coefficient corre- 
sponding to the m-th term in Eq. (4). Using the fact that the discretization scale a should 

be greater than the cutoff scale (1/k), wc have < P < 1. 

(b) The lattice is terminated at a large but finite A^. An intermediate point n is chosen, 
which we call the horizon TZ {— an), that separates the lattice points between the inside and 
outside. The horizon acts as a boundary of the bipartite system, (c) We confine our interest 
to the entanglement between the degrees of freedom inside and outside the horizon. The 

^ Such Hamiltonians are commonly encountered in liquid crystals and strain dynamics [20] by rewriting the 
above Hamiltonian in terms of the strain vector = Vj(^. 

^ Note that Kij contains nearest-neighbor (nn), next to nn (nnn), next to nnn (nnnn) and next to nnnn 
(nnnnn) interaction terms which originate in the derivative terms in Eq.(2). The central difference scheme 
we use is superior compared to the mid-point discretization scheme in the literature [4]. 
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reduced density matrix p^^^^^^^ (see Appendix A) is obtained by tracing over all the lattice 
points inside the horizon for the ground state wave-function of the discretized Hamiltonian 
(4). The resulting p^^^^^^^ is a mixed state of a bipartite system. Entanglement is computed 
as the von Neumann entropy associated with the reduced density matrix p^^^^^^^ [1, 2]: 

^ Tr reduced "^"^ reduced} 

In the next section, we show the results of numerical computations of entanglement entropy 
(following the method described in Appendix A) for our model. 

III. NUMERICAL RESULTS 

We compute the entanglement entropy for the discretized Hamiltonian (4) for different 
couphng strength P (Appendix B proves convergence of total entropy as i ^ oo). The 
computations are done using Matlab for the lattice size N — 300, 50 < n < 295 and the 
relative error in the computation of the entropy is 10~^. The computations are done for the 
following three different scenarios: 

(1) r = 0, e = 1: In Fig. (1), we have plotted InS*-^-* versus ln(7^/a) for different value 
of the coupling parameter P. The best fit (solid) curves in the two asymptotic regimes, 
P ^ and P ^ 1, show that the entropy scales approximately as area S ^ (TZ/a)"^. As we 
increase P, however, the prefactor increases. For P = 10~^ and 10~^, around the transition 
region between linear to non-linear, entropy increases by an order. 

(2) T = l,e = 0:In Fig. (2), \nS^^^ versus ln(7^/a) are plotted for different coupling 
parameter P, from which we infer the following: (i) For small P, as in the earlier case, 
entropy scales as area, (ii) As we increase P, above a critical value of P, say Pc (here 
Pc ~ 10~^) a new phase appears in the log — log plots of entropy vs n. (iii) For a given 
P > Pc, scahng of 5"^^^ changes drastically across some critical value of n, say ric (e.g. 
Tie ~ 139 for P = 10~^). For n > ric, S^^^ increases with n and scales approximately close 
to area. However, for n < tt-c, S'^^^ increases with decreasing n i.e. oc n " (a > 0)^ 
(iv) With increasing P, Uc increases. (V) Near P = Pc transition is sharp. For P > Pc, the 
transition region is wider. 

^ Note that as n — )■ 0, the entropy becomes zero. This is expected as the number of degrees of freedom 
gradually vanishes. 
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FIG. 1: Log-log plot of von Neumann entropy, S, versus the scaled radius of the sphere TZ/a = 
for r = 0, e = 1 and different values of coupling parameter P. The lattice size is = 300 and 
50 < n < 290. The dots represent the numerical output and solid lines denote lines of best fit. 
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FIG. 2: Log-log plot of von Neumann entropy, S, versus the scaled radius of the sphere IZ/a = n, 
for r = 1, e = and different values of coupling parameter P. The lattice size is = 300 and 
100 < n < 290. Uc for P = 10"^ 10"^, 10"^^ respectively, are 139, 167 and 197. 
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FIG. 3: Log-log plot of von Neumann entropy, S, versus the scaled radius of the sphere TZ/a = 
for r = 1, e = — 1 and different values of coupling parameter P. The lattice size is = 300 and 
100 < n < 290. The dots represent the numerical output and solid lines denote lines of best fit. 

(3) r = 1, e = —1: Fig. (3) also shows similar cross-over of entanglement entropy as in 
the earlier case. 

IV. UNDERSTANDING THE RESULTS 

The above numerical results clearly indicate that higher derivative terms seem to dras- 
tically modify the vacuum structure, due to which the entropy jumps by few orders of 
magnitude close to the critical point, and throws up a large number of interesting questions. 
In the rest of this article, we address three key questions related to this new phenomena: 

(I) What universal feature the higher derivative terms have on the scaling of the entropy? 

(II) What causes the entropy to increase by couple of orders of magnitude at the cross- 



over 



? 



(Ill) Does the sharp jump indicate phase transition? 
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FIG. 4: Log-Log plot of the entropy variation with only the 2nd order derivative term (left) and 
3rd order derivative term (right) present in the Hamiltonian. 

A. Inverse scaling 

To go about answering question (I), let us set the dispersion relation for all k modes to 
be co'^ — k^/hi^ and uJ^ — k^jhi^. Repeating the earlier analysis, In^*^^^ vs logn in Fig. (4) 
shows that the fourth order derivative term leads to the usual area law while the sixth order 
derivative term shows inverse scaling i.e. entropy increases as we decrease n. Comparing 
Fig. (4) with P — 1 plots in Fig. (2) show that increasing P leads to the dominance 
of the higher-order derivative terms implying that the inverse scaling of entropy uniquely 
corresponds to the presence of higher (sixth) order derivative terms. The system has two 
distinct entropy profiles, the area law due to (V0)^ and inverse scaling due to (V^0)^. Hence, 
for a particular coupling strength P^, the two distinct entropy profiles appear at two regimes 
and a crossover happens at a Uc that depends on P. 

So, why the sixth-order derivative term leads to inverse scaling while the fourth-order 
derivative does not? One could plausibly get a better understanding of the phenomena 
if we look at the density of states and the two-point correlation function (Appendix C) 
for the three dispersion relations as given in Table (I). The density of states for the two 
dispersion relations [k"^ and k'^) increase with energy while for uj"^ oc k^, it is constant. The 
density of states is a kinematical quantity, the two point-function (See Appendix C) contains 
information about the quantum fluctuations. The two point function is scale-invariant for 
the dispersion relation u"^ oc k^. We discuss the implications for condensed matter systems 
in conclusions. 
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Dispersion relation 


Density of states 


Two-point function 




D{uj)dLO oc io'^dio 




oc 


D{ijj)duj oc ijj^l'^du: 


G+ocl 


Lo^ oc k^ 


D{u})du} oc constant doj 


Scale- invariant 



TABLE I: Behavior of the density of states and two-point function for three dispersion relations: 
(i) bi / 0,62 = 63 = 0, (ii) 62 / 0,61 = 63 = 0, (iii) 63 / 0,61 = 62 = 0. 

B. Critical mode 

The plot of the distribution of entropy per partial wave {2i + l)^^ versus i can be used 
to answer question (II). Fig. (5) shows the entropy distribution for three values of n for 
e = 0, r = 1 and P = 10~^ — n = 200 (which falls in the positive slope region in Fig. 
2), n = 137 (which falls on the inverse scaling region in Fig. 2) and n = 110. It can be 
seen that, for n = 200, the contribution of large i falls off rapidly implying that only lower 
partial waves i contribute to the entropy. However, for n = 137, at a critical mode (ic), 
higher partial waves contribute significantly to entropy compared to lower i. For n = 110 
this happens at a lower i^- This implies that, for a fixed P, as n decreases, the turn around 
at which higher value of i contribution to the entropy occurs earlier which in-turn leads to 
inverse scaling of entropy. This feature can be derived analytically by looking at i ^ N 
limit of Hamiltonian (4) : 




FIG. 5: Log-Log plot of the entropy distribution per partial wave = {2i + l)Si and £ for n = 200, 
137 and 110 for P = 10"^ 
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(6) 



In the region when the higher-derivative term dominates, we have, 



(7) 



For a fixed P < 1, as seen in Fig. (5), as n increases the cross over occurs at larger i. In 
the case of P ~ 1, this imphes that the cross over occurs much earlier leading to inverse 
scaling of entropy^. For a detailed comparison of the entropy per mode among systems with 
different dispersion relations, see Appendix D. A qualitative understanding of excitation of 
higher modes (for P > Pc) can be understood using the quantum mechanical model of a 
particle in a box (see Appendix E). 

C. The cross-over 

1. Thermodynamic limit 

To answer question (III) that the change in the ground state signals phase transition, 
we first calculate the entanglement entropy for different values of N with P = 10~^,t = 
1, e = and accuracy 10~^ (see Appendix F). Fig. 6 clearly shows the ratio ric/N remains 
almost constant with increasing N, implying that, the cross-over is indeed a quantum phase 
transition. It is important to note there is a no phase transition for N — 75 indicating that 
Pc decreases with increasing N. 

2. Scaling transformations 

If Pc is the quantum critical point, from the Wilsonian view point, at this point the 
(quantum) fiuctuations occur at all length scales implying that action corresponding to 

Eq. (7) apparently implies that inverse scaling phase in entropy profiles should appear for any arbitrary 
low value for P > (n/ic)'^ . However, for very large £ the cntropic contribution from each £ is already 
suppressed and for low enough P, £c will become so large that contributions from £ > fc-modes will be 
too small to perturb the area scaling. 
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FIG. 6: Variation of vs iV for P = 10'^. 



Hamiltonian (2) or (4) is scale invariant. Under the scaling transformations (e 

= 1,2,3 (8) 

(9) 
(10) 

(11) 

the action corresponding to (4) is invariant for any value of the effective coupling strength 
P (since a and k scales inversely). Also, the action corresponding to (2) is invariant and k 
flows to zero. 
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V. DISCUSSIONS 



We have shown that linear, higher spatial derivative theories give rise to a new kind of 
quantum phase transition. The change in the scaling of entanglement entropy occurs due 
to two effects: 

1. Kinematical: Change in the density of states from increasing to constant. 

2. Dynamical: Excitation of larger i modes compared to the canonical theory. 

It is important to note that until now there is no fundamental understanding of the scaling 
of entanglement entropy with the density of states or the two-point correlation function. 
This work suggests a plausible link which will be explored in another work. 
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Recently, there has been interest in using entanglement as a probe to detect topological 
properties of many-body quantum states in two dimensional systems, like fractional quantum 
Hall effect [21, 22]. Recently, it is found that Graphene shows an unconventional sequence of 
fractional quantum Hall states [23]. Low-energy electronic states in Graphene are described 
by relativistic equation aka satisfying linear dispersion relation [24] . The higher-order deriva- 
tive term (62 7^ 0, 63 = 0) corresponds to Fermi gas has a constant density of state. Our 
work predicts that the above 2-D model will have an quantum phase transition. Gurrently, 
we are investigating the entanglement entropy to relate to the topological order. 

It is usually expected that the higher derivative terms tame the UV divergences. In the 
case of two coupled harmonic oscillators [4], it is known that the entropy diverges when the 
interaction strength is large. The infinitesimal boundary tends to increase the interaction 
strength leading to the divergence of the entanglement entropy [25]. As it is clear from 
our results, the divergence in the entanglement entropy can not be resolved by introducing 
higher derivative terms. 

Our results have interesting implications to black-hole entropy: Numerical results for 
semi-classical black-holes have shown that close to 90% of Hawking radiation is in s-wave 
[26] . Our results show that this is the case for the linear dispersion relation (or large black- 
hole limit). However large £ modes contribute, when the non- linear dispersion relations are 
dominant. It is also important to note that earlier analyzes of Trans-Planckian effects on 
the Hawking radiation have been performed for ^ = [16, 17]. In the light of our analysis 
those calculations have to be revisited. 

Due to Hawking radiation, the mass of the black-holes will decrease. As the size of the 
black-hole decreases, the curvature of the event-horizon increases and hence one need to 
include higher curvature terms [27-30]. Our analysis shows that as the size of the black-hole 
decreases, below a critical radius, the scaling of entropy changes instantaneously from area 
to inverse of area. It is interesting to look at the plausible implications of our result for the 
final stages of a microscopic black-hole. 

Our results also seem to be related to causal dynamical triangulation approach of 
Ambjorn et al [31]. There it has been shown that the continuum limit corresponds to the 
linear dispersion relation at large scales which in the UV limit changes to anisotropic scaling 
of space (with respect to time) corresponding to the dispersion relation F{k) oc k^. Recently, 
Horava [32] showed that gravity with dynamical critical exponent m = 3 is perturbatively 
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renormalizable (see also, [19]). Our analysis suggests that causal dynamical triangulation 
and Horava-Lifshitz model should also see a phase transition. These are under investigation. 
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Appendix 



Appendix A: Entanglement entropy from density matrix 



The Hamiltonian of a scalar field with modified dispersion relation, Eq. (2), in (3 + 1) 
dimensional flat space, is given by: 

1 



where k is the inverse length scale in the theory, e = 0, ±1 and r = 0, 1. 
Decomposing the field and its conjugate momentum in partial waves 



7r/m(r) 



Im 



Im 



yields: 



H = 



Im 



Im 

+ 



dr 



(Al) 



dr[ r r2 ^^'"^'^^ ^2 1 ^^2 



dV + 1) 



r |/<iV_lrfV_((l±i)*+3,(, + l)4V^'(' + l)^'P^ f 



\ dr^ r dr^ 



dr 



dr"^ 



(A2) 
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The discretizing scheme used in this work is different compared to that used by Srednicki [4]. 

Srednicki's mid-point discretization is not suited in the presence of higher-derivative terms. In this 

work, we have used central difference scheme, i.e., 

/(x + 1) - /(x - 1) 



fix) 
fix) 

which is 2nd order accurate. 



2Ax 

f{x + \)-2f{x) + f{x-\) 

(A.x)2 

J{x + 2) - 2J{x + 1) + 2f{x - 1) - fix - 2) 
2{Axf 



(A3) 



Eq. (A2) can be written as the Hamiltonian of a set of A'^ coupled harmonic oscillators 



1 ^ 

H = ^Hj = — 6ijTT] + ifj Kij ifi 

3 i,j 

where the off-diagonal elements of the matrix Kij represent the interactions: 
1 



(A4) 



K,. 
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^3 (2) + + 4 



^^^^^^ ^ a7JZ4 i^i'^^j'^ + '^i(iV-l)'^i(7V-l)] 



+ 
+ 

+ 

+ 



4A2 
(f + l) , 2e 



2z2 



+ 



+ 



Ff{N-l) + 



40,1 



(A -1)2 4 



a?K^ a^K^ \ 2 



^ + - 1) + + 1) + + 1) 



(z2- 1)2" 



^i3{i^l,N) 



2 + 



Z(Z + 1) 



+ 



4 4 



4 + ^(^ + 1) (4+ ^ 



U + 1)^ 



+ 



r (Fi{j-1) 



4 4 



(1 - + Fi{m{f)+ 



+ 



a2«;2 



4 + Z(Z + l)f^+ 1 



(i + l)2 



+ 



(j + 1) 
r rFi(i-l) 



(l-5a) + FiWF2W+ 



fc^F.W + + l)F,i^ + 1) + '^F.i^ + 1) - ^M!±^(l - .,.) 



(j + 1)^ r_ 

4j{j + 2) 02^2 a'^K 

(^ + 1) 



+ 1)^3 0" + 1) + 



(i + l)2--- ■ - 2 

Z(/ + l) , F2(j)-^2(i + 2) 



(J + 1)^ 



2 



r 



4i(i + 2) 02^2 a^K^ 
r 



4 4 



(i + 1)^ 

F3(z + l)-Fi(z + 2) 



^i,3+2 



Si,j-3 - ^((^i,j+4 + Si,j-(^f> 



where 



1 /(/ + 1) 2 3^0 + 1) 
Fi(z) = -1 - - - -^T^,F2{l) = ^ + ^3 ., 

- 1-^+ ^-2 ,i^4W--2 • 



(A6) 



13 



Note that Kij contains nearest-neighbor (nn), next to nn (nnn), next to nnn (nnnn) and next to 
nnnn (nnnnn) interaction terms due to the presence of higher derivative terms in (A2) . Schemati- 
cally, 

^ X X X X X ^ 

X X X X X X 
X X X X X X X 

xxxxxxxx 

xxxxxxxxx (A7) 
xxxxxxxx 

X X X X X X X 
X X X X X X 
X X X X X 

A brief description of how to calculate entropy from the above Hamiltonian is the following. The 
density matrix, tracing over the first n of N oscillators (r = r„+i, . . . ,rjv), is given by: 

/n 
Y[ <^(ri,...,r„;r„+i,...,rAr) </C*(ri,...,rn;r^+i,...,r^) (A8) 

i=l 

where r and r' represent radial distances of the points, outside the horizon, from the center. The 
ground state is 



V 



N 



1 



0(ri,...,rAr) = A^i exp(--fc|,.r,2). 



(A9) 



1=1 



(where T = Ur, U KU^ = Kd a diagonal matrix) the corresponding density matrix (A8) can be 
evaluated exactly (the superscript (0) signifies GS): 



Pieduced ~ exp [-(r^7r + r'^ir')/2 + r'^ (5 



where: 



A B \ 1 ^ 1 

B"^ C ' 2 



(AlO) 



(All) 



Note that B and /3 are non-zero if and only if there arc interactions. The Gaussian nature of the 

_i _i 

above density matrix lends itself to a series of diagonalisations [^7^"^ = diag, /3 = "yjy^VpV'^'j^^ , 
WPW^ = diag, Vi E V = W'^ {VjV'^)2VT], such that it reduces to a product of (N—n), 2-oscillator 
density matrices, in each of which one oscillator is traced over [4]: 



N-n 



n 



i=l 



vf + vf 



+ PiViv'i 



(A12) 
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The corresponding entropy is given by: 

N-n 



s=J2 (-Mi-^il- 34^in^^) 



(A13) 



Thus, for the full Hamiltonian H = Yliim ^im-, the entropy is: 

oo 

S = ^{21 + l)Si (A14) 

1=0 

where the degeneracy factor {21 + 1) follows from spherical symmetry of the Hamiltonian. In 

practice, we will replace the upper bound of the sum in the above to a large value Imax- For the 
interaction matrix without any correction i.e. with e = r = in (A5), the above entropy, computed 
numerically, turned out to be [3, 4]: 

5 = 0.3(n + 1/2)2 = 0.3 (^^j _ ^^^^-^ 

Appendix B: Asymptotic analysis: Convergence of entropic contribution 

For I » j, the Hamiltonian can be written as 

H ^ — 7r| + ^ I 1 + eP ^2 + ^ ) + interaction/perturbative terms... (Bl) 

In the case of e = 1 and r = in the Hamiltonian (Bl), we get ^{P ^ 0) = (where, 

] 4 4 

P = ^2^)) where as the corresponding result of Srednicki was [4] and ^(P — ^ 1) = jr- 

In the case of e = and r = 1, we get ^(P ^ 0) = ^ and ^(P ^ 1) = + ^ + These 
values confirm that the entropy converges for large I and justifies the use of an upper cutoff l^iax 
for numerical estimation. 



Appendix C: Calculation of correlation function 

The two point correlation function of the fields also provide information about structure of the 
scalar field from the linear to non-linear dispersion relations. The (equal time) two-point correlation 
function or the Wightman function is given by 

Sk exp[ifc.(x — y)] 

sin(fcr) dk (CI) 
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where r = \x — y\. Eq. (CI) can be rewritten for a general dispersion relation co = and with 
substitution kr = a as 

(C2) 



47r2 



a 



,m— 1 



Note that for linear (m = 1) and quadratic (m = 2) dispersion models correlation decays with 
increasing distance which explains the area-law behavior of the entropy. Interestingly when the 
third order correction (m = 3) is dominant, the correlation function essentially becomes scale 
invariant. 

Appendix D: Entropy per partial wave for different dispersion models 
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FIG. 7: Entropy distribution per partial wave FIG. 8: Entropy distribution per partial wave 
(2£ + 1)5^ and in linear dispersive system (21 + 1)S( and in quadratic dispersive sys- 
(m = 1), for n = 200 and n = 100. tem (m = 2), for n = 200 and n = 100. 

Fig. 7 and Fig. 8 shows that in linear and quadratic dispersive theories, with decreasing n, 
entropy per partial wave decreases. On the other hand for cubic dispersion theory exactly opposite 
behavior is found (see Fig. 9). 

Appendix E: Particle in a box model 



The Schrodinger equation corresponding to a general non-linear dispersion model is: 



dx^ 



(El) 



16 



n=200 



n=137 



n=100 
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FIG. 9: Entropy distribution per partial wave (2^ + 1)5^ and ^, in cubic dispersive system (m = 3), 
for n = 200, n = 137 and n = 100. 



The most general solution is given by 



2m 



Cj exp[(— l)^^£^2;^a;] , 



(E2) 



where Cj's are constants to be determined by the boundary conditions. In particular, for m = 1, 
2 and 3, the general solution is given by: 



m=l 



ci cos(£' 2 x) + C2 sin(i? 2 x) 



(E3) 



V'm=2 = exp 



£■40:; 



ci cos — -p^ + C2 sm 



^/2 



x/2 



+ exp 



V'm=3 = exp 



"71" 



(Six] . I Six 
C3 cos — -p^ + C4 sm 



v/2 



V2 



(E4) 



ci cos 



— — + C2 sm — — 



+ exp 



\/3£:i 



6X 



C3 cos 



<f ex\ . (Sex 

— — + C4 sm — — 



+C5 cos (^^ ^ + C6 sin 6 x^ 



(E5) 



For a particle in a box model, only exponentially growing and decaying solutions exist for even 



m, while for odd m one has stationary solutions [~ sin(f2mx)]. With the appropriate boundary 
conditions at x = and L, one finds the energy eigenvalues are Ei,: 

— K— J i/ = l,2,.. (E6) 

Ground state energy eigenvalue for m = 1 Lk "^^"^ 
Ground state energy eigenvalue for m > 1 tt 



which leads to: 
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-m=1 
m=3 




FIG. 10: Energy eigenvalues for m = 1 and 3, L = tt and varying k. 



For m> 1, with decreasing k, the ground state energy eigenvalue of the system increases. Fig. 
(10) shows that at k = ir/L, it crosses the ground state energy eigenvalue for the system satisfying 
linear dispersion relation (m = 1) and becomes more than that for k < tt/L. 

Eq. (E7) implies that, for k < n/L, the ground state energy eigenvalue in the non-linear 
dispersion theory is higher compared to that of linear dispersion theory. As in the field theory 
model, where there is a cross-over from the linear to non-linear regime, with increasing P (or 
decreasing k), for P > Pc, the system needs to readjust in such a way that the ground state energy 
of the system increases. In other words, the cross-over of the dispersion relation catalyzes larger 
population of higher energy quantum states compared to low-energy states. 

Appendix F: Transition in the thermodynamic limit 

Fig. 11 shows that, for fixed P with increasing N the transition happens approximately at 
the same ric/N though it looses it's sharpness. The inverse scaling does not appear for N = 75 
implying that Pc decreases with increasing N. 
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